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Thermoelectric response near a quantum critical point of /3-YbAlB4 and YbRh2Si2 
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The thermoelectric coefficients have been measured on the Yb-based heavy fermion compounds 
/3-YbAlB4 and YbRh2Si2 down to a very low temperature. We observe a striking difference in the 
behavior of the Seebeck coefficient, S in the vicinity of the Quantum Critical Point (QCP) in the two 
systems. As the critical field is approached, S/T enhances in /3-YbAlB4 but is drastically reduced in 
YbRh2Si2. While in the former system, the ratio of thermopower-to-specific heat remains constant, 
it drastically drops near the QCP in YbRh2Si2. In both systems, on the other hand, the Nernst 
coefficient shows a diverging behavior near the QCP. The results provide a new window to the way 
various energy scales of the system behave and eventually vanish near a QCP. 

PACS numbers: 71.27.+a, 72.15.Jf, 74.40.Kb 
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Heavy fermion compounds have been intensively stud- 
ied for their ability to use pressure or magnetic field to 
access a quantum critical point (QCP) from a magnetic 
phase or a Fermi liquid (FL) state [l|. One consequence 
of the QCP is an emergence of non-Fermi liquid (NFL) 
behavior revealing strong deviations from the FL. The 
behavior has been characterized by a set of exponents of 
the temperature dependence of the physical properties, 
such as the specific heat, the resistivity, and the suscep- 
tibility in the framework of the spin fluctuation theo- 
ries [2|-y|. So far, these theories have been succeeded to 
describe the NFL properties of a number of materials Q . 
Recently, however, there has been found a few hecwy 
fermion materials, e.g., CeCui-^Aua; Q, YbRh2Si2 01, 
and /3-YbAlB4 [1], exhibiting the NFL behavior in the 
neighborhood of the QCP for which the standard theory 
is insufficient. In spite of the considerable efforts, both 
experimental and theoretical, up to date nature of their 
unconventional criticality has still been debated. One in- 
gredient preventing a clarification of the QCPs is due to 
a difficulty of precise determination of exponents for the 
NFL behaviors. A nuclear contribution to the specific 
heat, for example, tends to mask an intrinsic behavior at 
low temperatures, and necessarily leads to ambiguity. 

Recently, new insight on the QCP via the thermoelec- 
tric coeflacients, e.g., the Seebeck and Nernst coefficients, 
has been examined j^, • An advantage of these quanti- 
ties is that only itinerant quasiparticle (QP) excitations 
are measured, so thus they would be complementary to 
the specific heat: a measure of both itinerant and lo- 
calized excitations. Moreover, in contrast to the specific 
heat, the thermoelectric response is not determined by 
the number of itinerant electrons present in the system 
and its magnitude is directly set by their relevant energy 
scale. In this paper, we report on a comparative study 
of the thermoelectric response near the QCP of the two 



representative Yb-based unconventional quantum criti- 
cal materials, /3-YbAlB4 and YbRh2Si2. We find that 
while the Nernst response behaves in a similar manner in 
the two systems, the Seebeck coefficient and its relation 
to the specific heat are quite distinct in the two cases, 
allowing discrimination of the nature of their quantum 
criticality. 

/3-YbAlB4 is the first Yb-based heavy fermion super- 
conductor (Tc = 80 mK) Q. The striking feature of 
this system is an emergence of the NFL behaviors un- 
der zero field at ambient pressure. This implies that 
the system is located near the QCP without the exter- 



nal tuning parameter [11|. Coexistence of valence fluc- 
tuation and the NFL behavior is another feature of this 
material [12[. YbRh2Si2 is a heavy fermion compound 



showing the NFL properties above the antiferromagnetic 
(AF) transition at Tn = 70 mK 0, [H. Significantly, 
Tn is found to be continuously suppressed with a small 
magnetic field (-Be = 0.66 T || c-axis) to zero accessing 
the field-induced QCP. An existence of additional energy 
scale T* discussed as the crossover line of the Kondo de- 
struction is a puzzling feature of YbRh2Si2 IJ]. A key 



characteristics of both systems, which goes beyond the 
standard theory, is a diverging enhancement of the ef- 
fective mass and the uniform susceptibility towards the 
QCP 0, H, Stimulated by these striking similarity, 
the theoretical efforts have been devoted to understand 
the unconventional quantum criticality of the two sys- 
tems in a unified way 15l-ll7l| , but no consensus has been 
yet reached. 

Single crystals of /3-YbAlB4 and YbRh2Si2 were pre- 
pared using the aluminum and indium flux, respec- 
tively [H, [3l ■ We used two crystals with different quali- 
ties for /3-YbAlB4 (the residual resistivity ratio {RRR) is 
130 and 270 for sample jj 1 and ji 2, respectively), and one 
crystal {RRR 90) for YbRh2Si2. The thermoelectric 
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coefficients were measured by employing a steady-state 
method in a dilution refrigerator. The heat current q 
was injected parallel to the a6-plane, and the magnetic 
field was applied parallel to the c-axis on the sample with 
a size of typically 0.5 x 0.1 x 0.005 mm'^ for /3-YbAlB4 and 
1.6x2.2x0.05 mm^ for YbRh2Si2, respectively. The ther- 
mal contacts with resistance of < 10 mfi at room tem- 
perature were made by using a spot welding technique. 
The same contacts were used to measure the resistivity 
by a standard four-contact method. 

Figure [I{a) shows the temperature dependence of the 
Seebeck coefficient divided by temperature —S{T)/T un- 
der zero field and 5 T for the two different samples of 
/3-YbAlB4 (sample jjl and p). What is remarkable is 
that the —S{T)/T curves for sample p and p converge 
at T < 0.15 K under each field, even though the resid- 
ual resistivity is different by a factor of 2.5 (the inset 
of Fig. ma)). This indicates that the Seebeck coeffi- 
cient takes an identical value at sufficiently low tempera- 
ture, irrespective of the sample quality. Using the Boltz- 
mann equation, the Seebeck coefficient is expressed as 
S/T = -(7rV3)(fc|/e)(91nCT(e)/ae),=,p (a: the electri- 
cal conductivity). In the simple case of a free electron gas, 
the term {d\na{e)/de)e=eF is approximated by the en- 
ergy derivative of the scattering time, (91nr(e)/9e)c=ep. 
Moreover, in the zcro-tcmpcrature limit, it can be further 
simplified as [l^ 

1 6 e n \ 3 J 

by expressing T(e) ~ Toe'', where N, ep, and n arc the 
density of states, the Fermi energy, and the carrier con- 
centration, respectively. Therefore, the Seebeck coeffi- 
cient is proportional to the density of states per electron 
as T — >■ 0. This fact gives rise to an intriguing aspect 
of S/T that it is interrelated with the electronic specific 
heat coefficient C /T = 70 via the dimensionless ratio 
q = {S/T^o)NAye ~ ±1 {Nav- the Avogadro number), 
since both quantities are proportional to the density of 
states. It has been shown that this semi-universality is 
held within logarithmic accuracy in a wide variety of met- 
als including the heavy fermions [20j . However, an insight 
into this ratio in the vicinity of the QCP remains an issue 
of debate. 

With the above considerations in mind, let us look 
at variations of S{T)/T under fields. Figure [T] shows 
the —S{T)/T vs T plots under several magnetic fields 
for (b) /?-YbAlB4 (sample p) and (c), (d) YbRhsSis, 
respectively. In both systems, the Seebeck coefficient 
takes a negative sign as observed in other Yb-based heavy 
fermion compounds [2^. First, we focus on the results of 
/3-YbAlB4. Under zero field, -S{T)/T exhibits a loga- 
rithmic increase down to 0.2 K and reaches a large value 
—S/T ~ 6 nV/K"^, being two order of magnitude larger 
than the values for simple metals such as Cu {S/T ^ 
-30 nV/K^) [2^. The large S/T was also found in other 




FIG. 1; (color online), (a) Temperature dependence of the 
Seebeck coefficient -S(T)/T for two crystals of /3-YbAlB4 
(sample [tl and (12) under zero field and 5 T. Inset: zero-field 
resistivity vs T plot for jjl and jJ2. (b) —S{T)/T curves under 
several fields for [12. The -S{T)/T curves for YbRli2Si2 (c) 
below and (d) above the critical field of 0.66 T, respectively. 



heavy-fermion compounds 0, [13, [l^l and attributed to 
a large effective mass. On further decreasing the tem- 
perature, —S{T)/T shows a steep increase followed by a 
sudden drop at 80 mK due to the superconducting transi- 
tion. By contrast, under the field of 25 mT which is larger 
than the upper critical field [2lj but very close to the 
QCP, the drop disappears and —S{T)/T is found to con- 
tinuously increase down to the lowest temperature. By 
increasing the field B > 0.5 T, the low-temperature in- 
crease is suppressed and —S{T)/T approaches constant, 
indicating that the system becomes the FL. 

Next, we turn to the results of YbRh2Si2. The 
—S{T)/T curves taken under the fields below and above 
the critical field {Be = 0.66 T) are shown in Fig.[TJc) and 
[TJd), respectively. Under zero field, —S{T)/T shows a 
logarithmic increase at high temperature T > 0.2 K, the 
same as /3-YbAlB4. However, upon cooling, —S{T)/T 
peaks around 0.1 K which is slightly above Tn, and then 
turns to decrease. By applying the fields up to 1.0 T, 
the peaks broaden and the decrease of —S{T)/T be- 
comes more remarkable. It should be emphasized that, 
in contrast to /3-YbAlB4 , we find a considerable decrease 
of —S{T)/T under the critical field, while the logarith- 
mic behavior is still observed at T > 0.2 K. Theoreti- 
cally, such logarithmic variation was predicted both in 
the spin-density-wave (SDW) [lH and the Kondo break- 
down scenarios [l^l- Although the experimental results 
agree with these predictions at high temperature, they 
contradict at low temperature, indicating that the loga- 
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FIG. 2: (color online). Contour plots of the Seebeck coeffi- 
cient -S/T for (a) /3-YbAlB4 and (b) YbRli2Si2 in the field- 
temperature plane. Tn, Tfl, and T* represent the Neel or- 
dering temperature, the crossover temperature to the Fermi 
liquid state, and the crossover line where the Hall coefficient 
shows a rapid change, respectively p^ . 



rithmic divergence of S/T does not necessarily capture 
the intrinsic quantum criticality. Rather, the contrasting 
behavior found at sufficiently low temperature would be 
crucial for understanding the QCP. As further increas- 
ing the fields, —S{T)/T turns to increase below 60 mK 
at 1.5 T and shows a continuous increase down to the 
lowest temperature at 2.0 T. Above 3 T, the saturation 
behavior eventually appears after taking broad maxima, 
which roughly coincide with the crossover temperature 
to the FL state . These observations are qualitatively 
consistent with the results of S. Hartmann [lO| in which 
the magnetic field is applied perpendicular to the c-axis. 

To clearly demonstrate the contrasting behavior of 
S/T near the QCP, we depict the contour plots of 
-S{T)/T in the B-T plane for (a) /3-YbAlB4 and (b) 
YbRh2Si2 in Fig. [21 respectively. One immediately no- 
tices that in the vicinity of the zero-field QCP {B^ = T), 
-S{T)/T for ^-YbAlB4 is very large, while the one for 
YbRli2Si2 is strongly suppressed around the field- induced 
QCP [Be = 0.66 T). In addition, we find two enhanced 
regions in the plot of YbRh2Si2; one is a zero-field re- 
gion above Tn and the other is a zero-temperature region 
around 2 T. Moreover, the T* line runs through between 
two enhanced regions. It seems that the evolution of 
—S/T towards the QCP due to increasing QP mass is 
unexpectedly suppressed around the QCP. An interpre- 
tation of this suppression will be discussed later. On the 
other hand, the fact that the large —S/T survives up to 
as high as 3 T, away from the QCP, implies an additional 
source for the enhancement. The field- induced FM fluc- 



tuations detected by the NMR measurements [2Jj could 
be a candidate. We also mention that the persistent evo- 
lution of —S/T in the low-field side of T* line, which 
amounts to as much as the one in the FL state, seems 
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FIG. 3: (color online). Field dependence of the Seebeck coef- 
ficient —S/T at 50 mK (left axis) and the electronic specific 
heat coefficient 70 (right axis) [3, @| for (a) /3-YbAlB4 and 
(b) YbRh2Si2, respectively, (c) The dimensionless ratio |g| 
vs B- Be plot for /3-YbAlB4 (left axis) and YbRh2Si2 (right 
axis). The dashed lines are guide to the eye. 



contradict with the Kondo destruction scenario predict 



26|. 



ing the presence of the small Fermi surface [25, 

Let us further discuss the contrasting behavior of S/T 
by examining the dimensionless ratio q = {S/T^Q)Np^e. 
Figure [3| displays the field dependence of —S{B)/T at 
50 mK and 70(B) 0, S, ^ for (a) /3-YbAlB4 and (b) 
YbRh2Si2, respectively. It is clearly seen that 70 di- 
verges as B ^ Be in both systems, indicating a diver- 
gence of the QP mass on approaching the QCP 0, @1 ■ As 
expected from Eq. (1), -S{B)/T for /3-YbAlB4 rapidly 
grows in parallel with 70 (-B). On the contrary, although 
-S{B)/T trucks 70(B) at high field in YbRh2Si2, it turns 
to decrease as B — > Sc. Consequently, \q\ appears to be 
field-independent \q\ ^ 4 in /3-YbAlB4, while the one 
of YbRh2Si2 shows a considerable decrease towards the 
QCP from a nearly constant value {\q\ ^ 1) at high field 
(Fig.[3]Jc)). Therefore, the semi- universality of q is held in 
/3-YbAlB4 even at the QCP, while it becomes invalid for 
YbRh2Si2 on approaching its QCP. This striking differ- 
ence of q seems to point to distinct nature of the quantum 
criticality in these systems. 

It has been pointed out that the semi-universality of 
q is satisfied when the whole Fermi surface (FS) equally 
contributes to S/T and 70 [1^. This is a case for the 
QCP with Q = fluctuations where the QP mass is en- 
hanced on the entire region of the FS ( the whole FS is 
"hot" ) , leading to the uniform enhancement of S/T and 
7o. However, the semi- universality is invalid for the case 
of AF (Q 7^ 0) QCP because the mass enhancement oc- 
curs on some parts of the FS, "hot regions" connected 



4 



> 



-(a) 


^-YbAlB - 


_ II 

1 0.5T 










^■"-v* 


< 1 < 



0.1 



T{K) 




T{K) 



FIG. 4: (color online) . Temperature dependence of the Nernst 
coefficient -v{T)/T for (a) /3-YbAlB4 and (b) YbRh2Si2 un- 
der several magnetic fields, respectively. The —u{T)/T curves 
for /3-YbAlB4 at B < 2.5 T are sifted for clarity. 



by the AF ordering wavevector Q. In this case, since 
the transport is dominated by the cold electrons, S/T is 
diminished due to a reduction of the hot electrons con- 
tribution, while 7o is enhanced because it always has a 
contribution from the whole FS. This leads to imbalance 
of S/T and 70, and results in the reduction of |g| from 
order of unity. In fact, a considerable decrease of [gj near 
the QCP was also found in CeCoIn5 0, for which the 
AF (SDW) scenario well accounts for the QCP. Accord- 
ing to this interpretation, the unchanged (|g| ^ 4) and 
the diminishing \q\ {\q\ < 1) towards the QCP suggest 
dominant Q = and finite Q fluctuations in /3-YbAlB4 
and YbRh2Si2, respectively. Apparently, this finding is 
at odds with the theoretical predictions which attempt 
to understand the quantum criticality of the two systems 
in a unified way [iq, ITtI . 

In contrast to the Secbeck coefficient, the Nernst co- 
efficient v/T is found to diverge in the vicinity of the 
QCP in both systems. Figure Sl^a) and HJb) represent 
—v{T)/T as a function of temperature for /3-YbAlB4 
and YbRh2Si2, respectively, under several fields. In (3- 
YbAlB4, we find a diverging increase of —v{T)/T^ reach- 
ing as much as 7 /iV/K^T at 0.5 T, while the enhance- 
ment is rapidly suppressed by the field of i? > 1 T. We 
note that —v[T)/T is very large, being three order of 
magnitude larger than that of simple metals [2^. Since 
the Nernst coefficient is expressed as the energy deriva- 
tive of the Hall angle, v/T cx {diaxiO-a/ de)^^^^, v/T is 
very small if tan^jj is only weakly dependent on energy, 
which is a good approximation in simple metals [29|. 
On the other hand, in the case of the energy-dependent 
tan^H, i^/T can be obtained from the following simple 

expression, valid as T — >■ [2^: ^ ~ '^tiere ^ 

represents the carrier mobility. Therefore, the divergence 
of —v{T)/T is interpreted as a consequence of vanishing 
of ep, defined as the characteristic energy scale of the FL. 
A similar diverging enhancement of —v{T)/T is also ob- 
served under the critical field {Be = 0.66 T) in YbRh2Si2 
(Fig.m^b)), although the behaviors of —v{T)/T is rather 
complicated including a sign change. By further apply- 



ing the fields, the diverging behavior is gradually sup- 
pressed and peaks are formed above 2 T (downward ar- 
rows). These peak temperatures correspond to TpL in 
Ref. [7[ , indicating that the FL is restored below the peak 
temperatures. The sigh change of —v{T)/T suggests an 
existence of two FSs with dominant carriers of opposite 
signs in accordance with the band calculations [30| , which 
may result in the relatively small magnitude of —v{T)/T 
compared with /3-YbAlB4 due to canceling of opposite 
sign contributions to the Nernst signal. 

In summary, we present the low-tcmperature ther- 
moelectric coefficient measurements of /3-YbAlB4 and 
YbRh2Si2. Both systems display the diverging enhance- 
ment of the Nernst coefficient on approaching their 
QCPs, indicative of the vanishing of the Fermi energy. 
On the other hand, the significant deviation from the 
semi-universal ratio q of the Secbeck coefficient and the 
electronic specific heat is found in YbRh2Si2 in the vicin- 
ity of the QCP, while the universality is sustained in j3- 
YbAlB4. This striking difference suggests distinct type 
of fiuctuations underlying the quantum criticality; Q = 
and finite Q fluctuations for /3-YbAlB4 and YbRh2Si2, 
respectively. Further experiments on other quantum crit- 
ical materials in terms of q would enable systematic clar- 
ification of the quantum criticality. 

This work is supported by Grants-in-Aids (Nos. 
23340099, 23740263, 21684019) for Scientific Research 
from the Japanese Society for the Promotion of Science, 
a Grant-in- Aid for Scientific Research on Innovative Ar- 
eas "Heavy Electrons" (Nos. 20102006, 23102706) of 
the Ministry of Education, Culture, Sports, Science, and 
Technology (MEXT), and a Grant- in- Aid for the Global 
COE Program from the MEXT through the Nanosciencc 
and Quantum Physics Project of the Tokyo Institute of 
Technology. 



[1] 

[2] 

[3: 

[6: 
[s: 

[9 

[10: 
[11 

[12 
[13 
[14 
[15 
[16 



P. Gegenwart, Q. Si, and F. Steglich, Nature Phys. 4, 
186 (2008). 

T. Moriya and T. Takimoto, J. Phys. Soc. Jpn. 64, 960 
(1995). 

J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 
A. J. Millis, Phys. Rev. B 7183, 48 (1993). 
W. Knafo, S. Raymond, P. Lejay, and J. Flouquet, Na- 
ture Phys. 5, 753 (2009). 

H. v. Lohneysen et al., Phys. Rev. Lett. 72, 3262 (1994). 

P. Gegenwart et al., Phys. Rev. Lett. 89, 056402 (2002). 

S. Nakatsuji et al., Nature Phys. 4, 603 (2008). 

K. Izawa et al., Phys. Rev. Lett. 99, 147005 (2007). 

S. Hartmann et al., Phys. Rev. Lett. 104, 096401 (2010). 

Y. Matsumoto et al., Science 21, 316 (2011). 

M. Okawa et al., Phys. Rev. Lett. 104, 247201 (2010). 

J. Custers et al.. Nature 423, 524 (2003). 

S. Paschen et al., Nature 432, 881 (2004). 

J. Custers et al., Phys. Rev. Lett. 186402, 104 (2010). 

T. Misawa, Y. Yamaji, and M. Imada, J. Phys. Soc. Jpn. 



5 



78, 084707 (2009). 
[17] S. Watanabe and K. Miyake, Phys. Rev. Lett. 105, 
186403 (2010). 

[18] R. T. Macaluso et al., Chem. Mater. 19, 1918 (2007). 

[19] G. Knebel et al., J. Phys. Soc. Jpn. 75, 114709 (2006). 

[20] K. Behnia, D. Jaccard, and J. Flouquet, J. Phys.: Con- 
dens. Matter 16, 5187 (2004). 

[21] K. Kuga et al., Phys. Rev. Lett 101, 137004 (2008). 

[22] I. Paul and G. Kotliar, Phys. Rev. B 64, 184414 (2001). 

[23] K.-S. Kim and C. Pepin, Phys. Rev. B 81, 205108 (2010). 

[24] K. Ishida et al., Phys. Rev. Lett. 89, 107202 (2002). 

[25] Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Nature 
413, 804 (2001). 



[26] P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. 
Phys.: Condens. Matter 13, R723 (2001). 

[27] The value of 70 at the QCP is taken from the temperature 
dependence of the electronic specific heat, Cbi{T)/T = 70 
at 50 mK [i,|l3]. In addition, for /3-YbAlB4, Cei{T)/T at 
50 mK is estimated from the extrapolation by assuming 
Cei(r)/roc -InT. 

[28] K. Miyake and H. Kohno, J. Phys. Soc. Jpn. 74, 254 
(2005). 

[29] K. Behnia, J. Phys.: Condens. Matter 21, 113101 (2009). 
[30] G. A. Wigger et al., Phys. Rev. B 76, 035106 (2007). 



